Evidence of breakdown of the spin symmetry in diluted 2D electron gases 
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Recent claims of an experimental demonstration of 
spontaneous spin polarisation in dilute electroiigases Q 
revived long standing theoretical discussions UQ- In 
two dimensions, the stabilisation of a ferromagnetic fluid 
might be hindered by the occurrence of the metal- 
insulator transition at low densities 0. To circumvent 
localisation in the two-dimensional electron gas (2DEG) 
we investigated the low populated second electron sub- 
band, where the disorder potential is mainly screened by 
the high density of the first subband. This letter reports 
on the breakdown of the spin symmetry in a 2DEG, re- 
vealed by the abrupt enhancement of the exchange and 
correlation terms of the Coulomb interaction, as deter- 
mined from the energies of the collective charge and spin 
excitations. Inelastic light scattering experiments and 
calculations within the time-dependent local spin-density 
approximation give strong evidence for the existence of 
a ferromagnetic ground state in the diluted regime. 

Recently we have shown experimentally as well as the- 
oretically that at low temperatures and zero magnetic 
field a 2D electron gas realized in a GaAs single quantum 
well (SQW) undergoes a first-order phase transition, as 
result of intersubband exchange-correlation terms when 
the second electron subband becomes populated |^ . The 
situation of two occupied subbands resembles that of dou- 
ble layer structures, a system with a rich ground state 
phase diagram, where a collapse of the normal metal- 
lic state towards a spontaneous magnetic order at very 
low densities was also predicted |^ Q . Evidence of this 
first order phase transition is the abrupt renormalisation 
of the second subband energy and its sudden popula- 
tion with electrons 8] as observed by photoluminescence 
(PL) spectra. The first-order character of the transition 
was confirmed by self-consistent density-functional cal- 
culations, including exchange interactions exactly [^ llClj| . 
The exact-exchange theory predicted also a spontaneous 
breakdown of the spin symmetry to take place for a low 
population of the second subband. In this situation a new 
phase was found in PL spectra at low temperatures and 
at zero magnetic field . The behaviour of this phase by 
varying density, temperature and external magnetic field 
speaks for the formation of spin-polarised domains with 
different in-plane magnetisation, as expected for planar 
ferromagnets to minimise the stray field. Moreover, last 
year possible ferromagnetism was claimed to be observed 
in transport experiments on a GaAs heterostructure with 



a single occupied subband 0. This evidence, though, 
is not conclusive due to the effects of localisation men- 
tioned above. In this work we investigate the elementary 
excitations of the 2DEG, observing a strong enhance- 
ment of the exchange-correlation contributions occurring 
at very low population of the second subband. Density- 
functional calculations within the time-dependent local 
spin-density approximation (TDLSDA) confirm the exis- 
tence of such spin polarised phase. 

The 2DEG forms in modulation-doped 
GaAs/Alo.ssGao.erAs single quantum wells of 20, 
25 and 30 nm well thickness grown by molecular-beam 
epitaxy. The 2DEG has a mean mobility of about 
8 X 10^ cm^ /Vs and a density such that only the lowest 
subband is occupied, with Fermi energies Ep between 
25 and 28 meV at 4.2 K. The electron density can 
be changed by applying a dc bias between the 2DEG 
and a metallic back contact. Photoluminescence and 
inelastic light scattering spectra were obtained using 
low power densities in the range of 0.5 to 5 W/cm^. 
Light scattering spectra were recorded in backscattering 
geometry using a line focus in order to avoid heating of 
the 2DEG at the required laser powers. 

The set of spectra in Fig. illustrates the evolution of 
the energy of the charge (CDE), spin (SDE) and single- 
particle excitations (SPE), associated with the intersub- 
band transitions 0^1 and 1-^2, when varying the den- 
sity of the electron gas. The shift of the SPEqi to lower 
energies and of the SPE12 to higher energies are con- 
sequences of the renormalisation of the second subband 
with increasing occupation. Furthermore, the change in 
the energy separation between excitations indicates a de- 
pendence of the many-body corrections on electron den- 
sity [Hill. 

As explained in the Appendix, within the formalism of 
TDLSDA !l3| and in the long wavelength limit (g « 
for backscattering geometry) , the energies of the charge 
and spin-density excitation Huj^ can be written in the 
following compact way 
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with (5nj.i-|_i = Tii — Hi^i being the difference in the sub- 
band densities and -Ei.i+i = — Ei, where Ei are the 
subband energies. In the paramagnetic (P) phase, the in- 
dex = CDE, SDE describes the uncoupled charge and 
spin density excitations. According to standard nota- 
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FIG. 1: Inelastic light scattering spectra of a 30 nm wide 
SQW obtained at 2 K, with different voltages and in polarised 
(II) and depolarised (_L) configuration. In the spectra with 
parallel polarisation the CDE associated with intersubband 
transitions between the first and the second subband — > 1, 
and between the second and the third subbands f — > 2 are 
apparent. The corresponding SDE, in contrast, are observed 
in depolarised configuration. The broader peak that appears 
in spectra with both polarisations is assigned to the single- 
particle excitation SPE and is centred at the intersubband 
energy Eij, with (ij) being (Of) or (12). 
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e{uj) is the contribution of the polar lattice to the di- 
electric function In. the ferromagnetic phase (F), 
charge and spin are intrinsically hnked, thus, there is a 
single collective mode of mixed character (denoted with 
subindex /i =M). The dcpolarisation shift a is a posi- 
tive Hartree contribution that in the paramagnetic phase 
typically dominates over the exchange-correlation term 
PcDE, whereas 7 is dominated by /3 in the ferromag- 
netic case. The explicit expressions for a, Pcde, Psde, 
and 7m in terms of matrix elements of the Hartree 
and exchange-correlation potentials are given in the Ap- 
pendix. 

In Fig. |21we show an unexpected behaviour of the spin- 
density mode: it becomes softer, i.e. decreases in abso- 
lute value when the Fermi level is resonant with the sec- 
ond subband, recovering the normal behaviour at higher 
densities. From the study of the formation of a non trivial 
magnetic order, intensively investigated in double-layer 
systems Is], a softening of the SDEqi can be consid- 
ered as an indication of ferromagnetism in the second 
subband. In our system the softening is rather weak, 
probably due to the large contribution of electrons in the 
paramagnetic first subband. Comparing the three QW 
widths (Fig.|2ta)) it can be seen that the SDE softening is 
enhanced for the narrowest QW. That speaks sensibly for 
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FIG. 2: The 0^1 spin-density mode softening as a function 
of the dimensionless parameter Ef/Eoi, being Ef/Eoi — 1 
the borderline between one and two occupied subbands. a) 
Comparison of the SDE softening among the 20, 25, and 30 
nm wide SQW. b) Comparison of the experimental SDE soft- 
ening (solid circles, left axis) for the 20 nm wide SQW with 
TDLSDA results (open circles, right axis). The lines are 
guides to the eye. 



a dependence of the possible ferromagnetic order on the 
strength of the wave function confinement. In Fig. Hfb) 
we compare the experimental data (solid points) with 
the excitation energies yielded by TDLSDA calculations 
(empty points). The left (right) axis corresponds to the 
experimental (calculated) results. The underestimation 
of the softening resulting from the calculations can be 
attributed to shortcomings of the TDLSDA in the treat- 
ment of exchange and correlation but also to the fact that 
the calculations always yield a paramagnetic ground state 
for this high total density. 

For this reason and in order to exclude the effect of the 
electrons in the first paramagnetic subband, we investi- 
gate higher intersubband excitations associated with the 
transition 1 — > 2. These excitations occur between the 
sparsely populated second subband and the empty third 
subband. Only electrons from the spin polarised phase 
in El contribute to the SDE12. In the series of spectra of 
Fig. nthe excitations SPE12, SDE12 and CDE12 become 
apparent when the second subband is populated. In par- 
ticular, the CDE12 shows up exactly on top of the SPE12 
due to the collapse of the Hartree term a^^ at very low 
densities In contrast, the SDE12 is always shifted 
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down from the SPE12. From the measured energies of 
the excitations obtained with cross polarisation and using 
Eq. 1^ we determine the many-body correction energy 
(27ii7^^). Multiplying by ni we avoid divergencies when 
the density goes to zero. The results are shown in Fig.O 
The striking result of this work is the observation of an 



> 



J 


■ 


2 K 


\ 

\ 





i\ 


\ 

\ 




10K 


■ ^ 

\ 

~ \ 




Vo 


\ 

\" 

\ 






\ 

V 






N 

■0 












- or 

1,1, 


1 


1 



12 3 

n, ( 10 cm ) 



FIG. 3: Enhancement of the exchange-correlation energy con- 
tribution 27117'^^, as a function of the density n\ of the sec- 
ond subband, for the 30 nm wide SQW and at temperatures 
T= 2K (circles), T= 5K (squares) and T= 10 K (triangles). 
Curves are a guide to the eye. 

enhancement of the exchange-correlation contribution at 
low temperature when the system enters the region of 
the spin-instability. This leads to a pronounced soft- 
ening of the collective mode associated with electronic 
transitions 1^2. The diluted electron gas of the second 
subband becomes thus unstable against spin-flip excita- 
tions, which trigger the transition into the ferromagnetic 
phase. By increasing the temperature a scaling down 
of the exchange-correlation energy occurs, in accordance 
with the disappearance of the new phase observed in the 
PL spectra, as discussed in previous work 

For a better understanding of the physical origin of this 
transition it is instructive to perform the corresponding 
calculation of the excitations. To avoid the effect of the 
high density in the first subband we constructed an aux- 
iliary structure, using for its intersubband energy Eq\ 
and density no, the experimental values E\i and ni, re- 
spectively. Doing so, the calculations can yield now a 
ferromagnetic ground state and the corresponding exci- 
tations. 

Figure 4 compares the experimental 7^^^ values, as ob- 
tained from the excitations measured in cross polarisa- 
tion at 2 K, with the calculated ones. TDLSDA predicts 
the occurrence of a ferromagnetic phase, characterised 
by 7m (open circles), for densities below 0.34x lO^^cm"^ 
and a paramagnetic phase otherwise, characterised by 
the excitonic shifts 13c de and Psde (triangles). Inspec- 
tion of the different terms in reveals that its strong 
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FIG. 4: Comparison between the calculated many-body cor- 
rections 7m (ferromagnetic solution), (3c de and Psde (para- 
magnetic solution), obtained within TDLSDA and yl^p from 
depolarised spectra at 2K of the 30 nm wide SQW. The 
dashed line divides the ferromagnetic (left) from the para- 
magnetic (right) region. 

increase at low densities is mainly due to the fact that 
the contribution from correlation has the opposite sign 
than the one to Psde in the paramagnetic phase. Thus, 
correlation adds to the effect of exchange in the spin- 
polarised case, which favours ferromagnetism. In fact, 
the calculated values of Psde and 7m at high and low 
densities, respectively, agree well with 7^^^. The range 
of stability of the ferromagnetic solution is also in nice 
agreement with the data of Fig. 3, where the strong 
enhancement of exchange-correlation vertex corrections, 
observed at low temperatures, occurs at densities be- 
low 1 X lO^-'^cm"^. In this way, we have obtained from 
light scattering experiments and local spin-density calcu- 
lations the first compelling evidence for the existence of 
a ferromagnetic ground state of 2DEGs in an effectively 
diluted regime; one of the exciting puzzles of many-body 
physics stated more than fifty years ago by F. Bloch 

APPENDIX A: METHODS 

The ground-state calculations have been performed 
in the local spin-density approximation (LSDA) in the 
framework of density-functional theory. We adjusted the 
geometry of the SQW in order to reproduce the Fermi 
energy Ep, and the two subbands energies Eq and Ei as 
obtained from photoluminescence experiments. 

The energy of the collective excitations of charge 
(CDE) and spin (SDE) were obtained with TDLSDA 
from the zeros of the secular determinant of the spin- 
dependent dielectric tensor (dete = 0), which can be 
written in terms of the Hartree and exchange-correlation 
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potentials in the subband representation The corre- 
sponding potential matrix elements, as a function of the 
wavevector q of the excitations, are 



(9) = 'aJ '^^J dz' (2Tre'^ / sq) exp{-q \z - z'\) 
X (z) (z)</)*^,(z') (^'), (Al) 



X [iC^j^cl^;) + (cr + cj')Jr,c{z) + (tct7^c(2;)] • 

(A2) 

(z) are the self-consistent wave functions that di- 
agonalise the effective one-dimensional LSDA Hamilto- 
nian, after assuming translational invariance along the 
{x,y) plane (area A). ,j' are subband indexes, and 

(J, a' = 1(— 1) for spin f (|). Eq. (|A1|) corresponds to 
the Hartree contribution, with e being the dielectric con- 
stant of GaAs. Eq. (|A2|) corresponds to the exchange- 
correlation contribution, with K, J, and / being deriva- 
tives of the exchange-correlation energy functional with 
respect to the density and magnetisation of the 2DEG ■ 



Making a two-subband approximation (z 



I ,J 



f), 



and in the long wavelength limit {q 0), the eigen- 
mode equation det e = can be solved analytically for 
{hwY . The resulting equations for the P phase take 
the form of Eq. J^l, provided the following identifica- 
tion is made: a = Vh,Pcde = — {V^c + KJi) /2, and 
PsDE = — (VjJ — V^Ji) /2. For the sake of simphcity, we 
have eliminated subband indexes, and we have used the 
symmetry Vh = V^^ = V^^ = Vjj^ = V^^ . In the F 
phase, 7m = {V^^ + VjJ). This single collective mode 
has a mixed CDE and SDE character, and accordingly it 
has both Hartree and exchange-correlation contributions, 
the latter being the dominant term in the F phase. 
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